This webpage will show the results of a bat survey study done in the Plumas National Forest in North California. The objective of this study is to determine the distribution of the different species of bats within the park. In order to do that we have performed occupancy models for the species present in the park. The results of this models will be shown as maps showing the probability of occurence of bats in each point, that is, if you see a value of 1, there is a 100% chance of finding a bat in that point, if there is a value of 0 there is 0% chance of finding that specie in that point, if there is a value of 0.5 there is a 50% chance of finding that specie in that point.
Another result
In this area 0 means absence, and 1 means prescence. This table has for each site (ID), every specie and day, so for example if Mylu1=0, that means that for Myotis lucifugus (common name Little Brown bat, was detected on day one for that particular site).
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | Model 7 | Model 8 | |
|---|---|---|---|---|---|---|---|---|
| psi(Int) | -1.23 | -1.22 | -1.32 | -1.41 | -1.68 | -1.80 | -1.23 | -1.59 |
| (1.04) | (1.06) | (1.06) | (1.05) | (1.10) | (1.12) | (1.04) | (1.09) | |
| psi(Burn.intensity.basal) | 151.57 | 145.77 | 131.30 | |||||
| (203.26) | (200.99) | |||||||
| p(Int) | -3.84*** | -4.15*** | -3.93*** | -3.48* | -3.89*** | -3.92*** | -3.23* | -3.75** |
| (1.17) | (1.23) | (1.19) | (1.35) | (1.18) | (1.18) | (1.32) | (1.41) | |
| p(Meantemp) | 0.13 | 0.15* | 0.14 | 0.16* | 0.14 | 0.14 | 0.15 | 0.19* |
| (0.07) | (0.07) | (0.07) | (0.08) | (0.07) | (0.07) | (0.08) | (0.08) | |
| psi(Burn.intensity.Canopy) | 102.87 | 185.56 | 152.69 | |||||
| (199.23) | (208.80) | (16777.81) | ||||||
| psi(Burn.intensity.soil) | 69.41 | 166.50 | 163.31 | 116.72 | ||||
| (205.06) | (16777.82) | (198.86) | ||||||
| p(sdtemp) | -0.29 | -0.31 | -0.34 | |||||
| (0.31) | (0.31) | (0.33) | ||||||
| Log Likelihood | -33.43 | -33.65 | -33.87 | -32.74 | -32.89 | -32.91 | -33.08 | -33.10 |
| AICc | 75.83 | 76.28 | 76.73 | 76.99 | 77.29 | 77.32 | 77.66 | 77.70 |
| Delta | 0.00 | 0.46 | 0.90 | 1.16 | 1.46 | 1.49 | 1.84 | 1.87 |
| Weight | 0.07 | 0.05 | 0.04 | 0.04 | 0.03 | 0.03 | 0.03 | 0.03 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||||||||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | Model 7 | Model 8 | Model 9 | Model 10 | Model 11 | Model 12 | Model 13 | Model 14 | Model 15 | Model 16 | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| psi(Int) | 1.02** | 0.58 | 0.57 | 1.02** | 0.60 | 0.54 | 1.06** | 0.67 | 0.65 | 0.67 | 0.67 | 0.70 | 0.49 | 1.03** | 0.58 | 0.58 |
| (0.35) | (0.44) | (0.43) | (0.35) | (0.44) | (0.45) | (0.37) | (0.43) | (0.43) | (0.43) | (0.44) | (0.44) | (0.44) | (0.36) | (0.44) | (0.44) | |
| p(Int) | 2.29** | 2.30** | 3.12** | 3.07** | -0.88 | -1.01 | -0.73 | 2.30** | 3.12** | 3.12** | -0.87 | -0.87 | 0.43 | 0.98 | 0.76 | -0.46 |
| (0.77) | (0.77) | (0.98) | (0.98) | (0.82) | (0.79) | (0.81) | (0.77) | (0.98) | (0.98) | (0.82) | (0.82) | (1.47) | (1.61) | (1.63) | (0.90) | |
| p(Meanhum) | -0.02* | -0.02* | -0.03* | -0.02* | -0.02* | -0.03* | -0.03* | -0.01 | -0.02 | -0.01 | ||||||
| (0.01) | (0.01) | (0.01) | (0.01) | (0.01) | (0.01) | (0.01) | (0.01) | (0.01) | (0.01) | |||||||
| psi(Burn.intensity.soil) | 0.37 | 0.41 | 0.44 | 8.55 | 7.31 | 0.41 | 0.48 | |||||||||
| (0.27) | (0.29) | (0.32) | (6.13) | (4.85) | (0.30) | (0.33) | ||||||||||
| p(sdhum) | -0.06 | -0.06 | -0.06 | -0.06 | -0.04 | |||||||||||
| (0.04) | (0.04) | (0.04) | (0.04) | (0.04) | ||||||||||||
| p(Meantemp) | 0.12* | 0.12* | 0.11 | 0.12* | 0.12* | 0.08 | 0.06 | 0.07 | 0.11* | |||||||
| (0.06) | (0.05) | (0.06) | (0.06) | (0.06) | (0.06) | (0.07) | (0.07) | (0.06) | ||||||||
| psi(Burn.intensity.Canopy) | -6.74 | 0.28 | 0.34 | 0.36 | ||||||||||||
| (4.96) | (0.25) | (0.28) | (0.31) | |||||||||||||
| psi(Burn.intensity.basal) | 0.24 | 0.25 | -4.20 | |||||||||||||
| (0.21) | (0.23) | (2.84) | ||||||||||||||
| Log Likelihood | -82.47 | -81.35 | -80.11 | -81.39 | -81.53 | -80.29 | -82.86 | -81.67 | -80.41 | -80.50 | -81.87 | -81.96 | -79.41 | -82.06 | -80.80 | -80.87 |
| AICc | 171.52 | 171.68 | 171.72 | 171.76 | 172.03 | 172.09 | 172.28 | 172.32 | 172.32 | 172.49 | 172.72 | 172.90 | 172.98 | 173.09 | 173.11 | 173.25 |
| Delta | 0.00 | 0.16 | 0.20 | 0.24 | 0.51 | 0.57 | 0.77 | 0.81 | 0.81 | 0.98 | 1.20 | 1.38 | 1.46 | 1.58 | 1.59 | 1.73 |
| Weight | 0.03 | 0.03 | 0.03 | 0.03 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.01 | 0.01 | 0.01 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||||||||||||||||
| Model 1 | Model 2 | Model 3 | |
|---|---|---|---|
| psi(Int) | -2.23* | -2.41* | -2.34* |
| (0.98) | (1.04) | (1.02) | |
| psi(Burn.intensity.basal) | -75.21 | -138.72 | -122.96 |
| (68.15) | (98.83) | (91.89) | |
| psi(Burn.intensity.Canopy) | 121.93 | 226.24 | 200.39 |
| (111.94) | (162.42) | (151.00) | |
| psi(Burn.intensity.soil) | -20.25 | -39.44 | -34.72 |
| (20.71) | (30.09) | (27.97) | |
| p(Int) | 3.56* | 4.34* | 3.60* |
| (1.44) | (1.86) | (1.57) | |
| p(Meantemp) | -0.31** | -0.33** | -0.31** |
| (0.10) | (0.11) | (0.10) | |
| p(sdhum) | -0.05 | ||
| (0.08) | |||
| p(sdtemp) | -0.01 | ||
| (0.28) | |||
| Log Likelihood | -31.97 | -30.93 | -31.29 |
| AICc | 78.08 | 78.80 | 79.53 |
| Delta | 0.00 | 0.72 | 1.45 |
| Weight | 0.17 | 0.12 | 0.08 |
| Num. obs. | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | |||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | |
|---|---|---|---|---|---|
| psi(Int) | -0.04 | -0.45 | -0.73 | 0.38 | 0.05 |
| (0.95) | (0.71) | (0.61) | (1.15) | (1.01) | |
| p(Int) | 4.10 | -0.28 | -2.14 | 2.53 | 4.92 |
| (4.32) | (1.06) | (1.70) | (3.89) | (4.28) | |
| p(Julian) | -0.03 | -0.02 | -0.03 | ||
| (0.02) | (0.02) | (0.02) | |||
| p(Meanhum) | -0.01 | -0.01 | |||
| (0.01) | (0.01) | ||||
| p(Meantemp) | 0.09 | ||||
| (0.14) | |||||
| psi(Burn.intensity.basal) | -33.38 | ||||
| (30.21) | |||||
| psi(Burn.intensity.Canopy) | 44.81 | ||||
| (40.51) | |||||
| Log Likelihood | -40.67 | -40.92 | -41.12 | -38.70 | -40.21 |
| AICc | 87.92 | 88.42 | 88.81 | 88.90 | 89.40 |
| Delta | 0.00 | 0.50 | 0.89 | 0.98 | 1.47 |
| Weight | 0.11 | 0.08 | 0.07 | 0.07 | 0.05 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | |||||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | Model 7 | Model 8 | Model 9 | Model 10 | Model 11 | Model 12 | Model 13 | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| psi(Int) | -0.93* | -0.92* | -0.95* | -0.45 | -0.95* | -0.44 | -0.53 | -0.84 | -0.93* | -0.85 | -0.81 | -0.94* | -0.53 |
| (0.45) | (0.44) | (0.46) | (0.34) | (0.44) | (0.35) | (0.32) | (0.45) | (0.43) | (0.46) | (0.45) | (0.45) | (0.33) | |
| psi(Burn.intensity.Canopy) | 0.42 | 0.31 | 0.31 | ||||||||||
| (0.25) | (0.20) | (0.22) | |||||||||||
| p(Int) | 6.21 | 6.33 | 6.95 | 5.25 | -3.53* | 7.10 | -3.14 | 8.20 | -3.57* | 8.62 | 8.26 | -3.47* | -1.18 |
| (5.19) | (5.19) | (5.27) | (5.27) | (1.69) | (5.25) | (1.61) | (5.38) | (1.70) | (5.54) | (5.40) | (1.67) | (1.10) | |
| p(Julian) | -0.06* | -0.06* | -0.06* | -0.05 | -0.05 | -0.05 | -0.06 | -0.05 | |||||
| (0.03) | (0.03) | (0.03) | (0.03) | (0.03) | (0.03) | (0.03) | (0.03) | ||||||
| p(Meantemp) | 0.22* | 0.22* | 0.22* | 0.19 | 0.20 | 0.18 | 0.20 | 0.19 | |||||
| (0.11) | (0.11) | (0.11) | (0.11) | (0.11) | (0.11) | (0.11) | (0.11) | ||||||
| p(sdtemp) | 0.57 | 0.59 | 0.58 | 0.56 | 0.56* | 0.73 | 0.56 | 0.74* | 0.57* | 0.75* | 0.75* | 0.57* | 0.68 |
| (0.32) | (0.32) | (0.32) | (0.33) | (0.27) | (0.37) | (0.29) | (0.36) | (0.27) | (0.36) | (0.37) | (0.28) | (0.35) | |
| psi(Burn.intensity.basal) | 0.32 | 0.22 | 0.22 | ||||||||||
| (0.20) | (0.15) | (0.17) | |||||||||||
| psi(Burn.intensity.soil) | 0.48 | 0.35 | 0.31 | ||||||||||
| (0.29) | (0.26) | (0.22) | |||||||||||
| Log Likelihood | -48.87 | -48.89 | -48.96 | -50.73 | -50.88 | -52.15 | -52.15 | -50.95 | -50.96 | -51.02 | -51.04 | -51.09 | -53.58 |
| AICc | 111.90 | 111.94 | 112.08 | 112.97 | 113.26 | 113.27 | 113.28 | 113.39 | 113.42 | 113.53 | 113.59 | 113.68 | 113.72 |
| Delta | 0.00 | 0.04 | 0.18 | 1.06 | 1.35 | 1.37 | 1.38 | 1.49 | 1.52 | 1.63 | 1.69 | 1.78 | 1.82 |
| Weight | 0.04 | 0.04 | 0.04 | 0.03 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | |||||||||||||
| Model 1 | Model 2 | Model 3 | |
|---|---|---|---|
| psi(Int) | 3.62 | 4.11 | 16.05 |
| (12.36) | (18.09) | (17.95) | |
| p(Int) | 7.34 | 6.54 | 6.82 |
| (4.13) | (4.30) | (4.22) | |
| p(Julian) | -0.05* | -0.06* | -0.05* |
| (0.02) | (0.02) | (0.02) | |
| p(Meantemp) | 0.08 | ||
| (0.08) | |||
| psi(Burn.intensity.Canopy) | 18.92 | ||
| (23.46) | |||
| psi(Burn.intensity.soil) | -26.43 | ||
| (32.33) | |||
| Log Likelihood | -35.47 | -35.01 | -33.82 |
| AICc | 77.51 | 79.00 | 79.13 |
| Delta | 0.00 | 1.49 | 1.62 |
| Weight | 0.08 | 0.04 | 0.03 |
| Num. obs. | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | |||
| Model 1 | Model 2 | Model 3 | Model 4 | |
|---|---|---|---|---|
| psi(Int) | 0.68 | 0.71 | 0.75 | 1.31** |
| (0.51) | (0.49) | (0.54) | (0.45) | |
| psi(Burn.intensity.basal) | -8.76 | -28.35 | ||
| (10.44) | (52.78) | |||
| psi(Burn.intensity.Canopy) | 13.29 | 44.74 | ||
| (16.49) | (85.78) | |||
| p(Int) | 1.22* | 1.25* | 2.60* | 1.13 |
| (0.61) | (0.62) | (1.24) | (0.63) | |
| p(Meanhum) | -0.01 | -0.01 | -0.02* | -0.01 |
| (0.01) | (0.01) | (0.01) | (0.01) | |
| psi(Burn.intensity.soil) | 0.70 | |||
| (0.56) | ||||
| p(Meantemp) | -0.07 | |||
| (0.05) | ||||
| Log Likelihood | -85.66 | -87.18 | -84.87 | -89.00 |
| AICc | 182.81 | 183.33 | 183.90 | 184.57 |
| Delta | 0.00 | 0.52 | 1.09 | 1.76 |
| Weight | 0.13 | 0.10 | 0.07 | 0.05 |
| Num. obs. | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||||
| Model 1 | Model 2 | Model 3 | |
|---|---|---|---|
| psi(Int) | -0.83 | -1.14 | -1.05 |
| (0.84) | (0.74) | (0.78) | |
| psi(Burn.intensity.Canopy) | 72.60 | 18.52 | 30.46 |
| (130.37) | (13.35) | (40.77) | |
| psi(Burn.intensity.soil) | -86.32 | -21.35 | -35.65 |
| (156.22) | (15.76) | (48.94) | |
| p(Int) | 13.69* | 7.91 | 8.69 |
| (6.07) | (4.67) | (5.05) | |
| p(Julian) | -0.07* | -0.05 | -0.06* |
| (0.03) | (0.03) | (0.03) | |
| p(Meantemp) | -0.20 | ||
| (0.12) | |||
| p(Meanhum) | 0.02 | ||
| (0.02) | |||
| Log Likelihood | -25.76 | -27.54 | -26.35 |
| AICc | 65.68 | 66.57 | 66.85 |
| Delta | 0.00 | 0.90 | 1.18 |
| Weight | 0.19 | 0.12 | 0.10 |
| Num. obs. | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | |||
| Model 1 | Model 2 | |
|---|---|---|
| psi(Int) | -1.93* | -1.99* |
| (0.79) | (0.79) | |
| psi(Burn.intensity.basal) | 42.88 | 50.65 |
| (29.30) | (35.10) | |
| psi(Burn.intensity.Canopy) | -78.13 | -92.09 |
| (51.58) | (62.19) | |
| psi(Burn.intensity.soil) | 24.28 | 28.39 |
| (14.59) | (17.88) | |
| p(Int) | -5.56** | -6.56** |
| (1.76) | (2.20) | |
| p(Meantemp) | 0.32* | 0.35** |
| (0.13) | (0.13) | |
| p(sdtemp) | 0.22 | |
| (0.27) | ||
| Log Likelihood | -26.59 | -26.17 |
| AICc | 67.34 | 69.29 |
| Delta | 0.00 | 1.95 |
| Weight | 0.14 | 0.05 |
| Num. obs. | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | Model 7 | Model 8 | Model 9 | Model 10 | Model 11 | Model 12 | Model 13 | Model 14 | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| psi(Int) | -0.92 | -0.27 | 0.27 | -0.85 | -1.31 | -1.13 | 0.06 | -1.08 | -1.11 | -1.63 | -0.73 | -0.65 | -0.69 | 0.23 |
| (1.35) | (1.61) | (2.14) | (1.07) | (1.16) | (1.46) | (1.88) | (1.46) | (1.07) | (1.24) | (1.36) | (1.37) | (1.40) | (1.84) | |
| psi(Burn.intensity.basal) | -51.10 | -51.32 | -23.85 | -18.75 | 0.78 | |||||||||
| (48.85) | (41.90) | (31.11) | (20.57) | (1.08) | ||||||||||
| psi(Burn.intensity.Canopy) | 70.07 | 70.52 | -34.60 | 0.91 | ||||||||||
| (66.72) | (57.35) | (37.55) | (1.11) | |||||||||||
| p(Int) | -4.73** | -4.38* | -8.63* | 6.23 | 5.17 | -4.39* | 0.46 | -4.26* | -2.19 | 5.54 | -9.13** | -9.15** | -9.10** | -0.07 |
| (1.77) | (2.13) | (3.76) | (7.28) | (7.13) | (1.74) | (6.52) | (1.73) | (1.76) | (5.77) | (3.54) | (3.50) | (3.52) | (7.79) | |
| p(Meanhum) | 0.03 | 0.03 | 0.05 | 0.03 | 0.03 | 0.03 | 0.03 | 0.05 | 0.05 | 0.05 | 0.05 | |||
| (0.02) | (0.02) | (0.03) | (0.02) | (0.02) | (0.02) | (0.02) | (0.03) | (0.03) | (0.03) | (0.03) | ||||
| p(Meantemp) | 0.16 | 0.05 | 0.19 | 0.19 | 0.19 | 0.17 | ||||||||
| (0.11) | (0.10) | (0.11) | (0.11) | (0.11) | (0.12) | |||||||||
| p(Julian) | -0.04 | -0.05 | -0.03 | -0.04 | -0.05 | |||||||||
| (0.04) | (0.04) | (0.03) | (0.03) | (0.04) | ||||||||||
| psi(Burn.intensity.soil) | 40.49 | 42.88 | 31.99 | 0.92 | ||||||||||
| (51.88) | (46.24) | (34.12) | (1.30) | |||||||||||
| Log Likelihood | -20.41 | -22.90 | -21.91 | -23.18 | -19.48 | -21.05 | -22.43 | -21.20 | -23.70 | -21.25 | -21.31 | -21.36 | -21.36 | -21.39 |
| AICc | 52.31 | 52.37 | 52.80 | 52.93 | 53.11 | 53.59 | 53.84 | 53.90 | 53.98 | 54.01 | 54.12 | 54.22 | 54.23 | 54.27 |
| Delta | 0.00 | 0.06 | 0.48 | 0.62 | 0.80 | 1.28 | 1.53 | 1.59 | 1.67 | 1.69 | 1.81 | 1.90 | 1.91 | 1.96 |
| Weight | 0.07 | 0.06 | 0.05 | 0.05 | 0.04 | 0.03 | 0.03 | 0.03 | 0.03 | 0.03 | 0.03 | 0.03 | 0.03 | 0.02 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||||||||||||||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | Model 7 | Model 8 | |
|---|---|---|---|---|---|---|---|---|
| psi(Int) | -1.07* | -0.97* | -0.95* | -0.95 | -0.86 | -1.08* | -1.02* | -0.83 |
| (0.49) | (0.46) | (0.46) | (0.52) | (0.50) | (0.49) | (0.50) | (0.49) | |
| psi(Burn.intensity.soil) | 0.53* | 0.51* | 0.53* | 0.53* | ||||
| (0.24) | (0.24) | (0.23) | (0.24) | |||||
| p(Int) | 2.02* | 2.03* | 2.04* | 3.94 | 4.02* | 1.61 | -2.14 | 3.95 |
| (1.01) | (1.01) | (1.01) | (2.01) | (2.01) | (1.07) | (4.77) | (2.02) | |
| p(Meanhum) | -0.03* | -0.03* | -0.03* | -0.04* | -0.04* | -0.03* | -0.03* | -0.04* |
| (0.01) | (0.01) | (0.01) | (0.02) | (0.02) | (0.01) | (0.01) | (0.02) | |
| psi(Burn.intensity.Canopy) | 0.45* | 0.44* | ||||||
| (0.21) | (0.22) | |||||||
| psi(Burn.intensity.basal) | 0.33* | 0.32 | ||||||
| (0.16) | (0.16) | |||||||
| p(Meantemp) | -0.11 | -0.12 | -0.11 | |||||
| (0.10) | (0.10) | (0.10) | ||||||
| p(sdhum) | 0.05 | |||||||
| (0.06) | ||||||||
| p(Julian) | 0.02 | |||||||
| (0.02) | ||||||||
| Log Likelihood | -57.63 | -57.84 | -57.91 | -57.02 | -57.17 | -57.18 | -57.24 | -57.30 |
| AICc | 124.25 | 124.66 | 124.79 | 125.54 | 125.84 | 125.87 | 125.98 | 126.10 |
| Delta | 0.00 | 0.41 | 0.55 | 1.30 | 1.59 | 1.62 | 1.74 | 1.85 |
| Weight | 0.06 | 0.05 | 0.05 | 0.03 | 0.03 | 0.03 | 0.03 | 0.02 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||||||||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | Model 7 | Model 8 | Model 9 | Model 10 | Model 11 | Model 12 | Model 13 | Model 14 | Model 15 | Model 16 | Model 17 | Model 18 | Model 19 | Model 20 | Model 21 | Model 22 | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| psi(Int) | -0.29 | -0.30 | -0.34 | -0.05 | -0.23 | -0.27 | -0.22 | -0.30 | -0.31 | -0.35 | -0.28 | -0.32 | -0.27 | -0.28 | -0.33 | -0.29 | 0.28 | 0.19 | 0.20 | 0.16 | 0.19 | -0.25 |
| (0.41) | (0.41) | (0.42) | (0.56) | (0.44) | (0.45) | (0.43) | (0.40) | (0.41) | (0.42) | (0.41) | (0.43) | (0.41) | (0.41) | (0.42) | (0.41) | (0.37) | (0.32) | (0.32) | (0.31) | (0.32) | (0.43) | |
| psi(Burn.intensity.basal) | 0.27 | -72.83 | 0.25 | 0.25 | 0.25 | 0.26 | ||||||||||||||||
| (0.16) | (85.88) | (0.15) | (0.15) | (0.15) | (0.15) | |||||||||||||||||
| p(Int) | 0.91 | 0.91 | 0.91 | 6.30* | 1.31 | 1.28 | 1.25 | 1.02 | 1.02 | 1.02 | 0.19 | 0.20 | 0.19 | 0.95 | 0.93 | 0.91 | 1.45 | 0.89 | 0.14 | 1.04 | 0.74 | 1.65 |
| (0.48) | (0.48) | (0.48) | (3.07) | (1.04) | (1.04) | (1.04) | (0.79) | (0.79) | (0.79) | (0.72) | (0.72) | (0.72) | (3.12) | (3.12) | (3.11) | (1.04) | (0.49) | (0.74) | (0.79) | (3.17) | (1.11) | |
| p(sdhum) | -0.04 | -0.04 | -0.04 | -0.03 | -0.04 | |||||||||||||||||
| (0.04) | (0.04) | (0.04) | (0.04) | (0.04) | ||||||||||||||||||
| psi(Burn.intensity.Canopy) | 0.36 | 102.16 | 0.35 | 0.34 | 0.34 | 0.35 | 0.36 | |||||||||||||||
| (0.21) | (120.18) | (0.21) | (0.20) | (0.21) | (0.21) | (0.22) | ||||||||||||||||
| psi(Burn.intensity.soil) | 0.39 | 0.38 | 0.38 | 0.38 | 0.39 | |||||||||||||||||
| (0.22) | (0.23) | (0.22) | (0.22) | (0.22) | ||||||||||||||||||
| p(Julian) | -0.02 | -0.00 | -0.00 | -0.00 | -0.00 | |||||||||||||||||
| (0.02) | (0.02) | (0.02) | (0.02) | (0.02) | ||||||||||||||||||
| p(Meantemp) | -0.15** | -0.06 | -0.06 | -0.06 | -0.07 | -0.06 | ||||||||||||||||
| (0.06) | (0.08) | (0.08) | (0.08) | (0.08) | (0.08) | |||||||||||||||||
| p(sdtemp) | -0.11 | -0.11 | -0.11 | -0.12 | ||||||||||||||||||
| (0.18) | (0.18) | (0.18) | (0.18) | |||||||||||||||||||
| p(Meanhum) | 0.01 | 0.01 | 0.01 | 0.01 | ||||||||||||||||||
| (0.01) | (0.01) | (0.01) | (0.01) | |||||||||||||||||||
| Log Likelihood | -73.03 | -73.04 | -73.06 | -70.50 | -73.18 | -73.20 | -73.22 | -73.27 | -73.27 | -73.28 | -73.29 | -73.30 | -73.30 | -73.44 | -73.44 | -73.45 | -74.78 | -74.82 | -74.97 | -74.97 | -75.16 | -72.76 |
| AICc | 155.03 | 155.05 | 155.10 | 155.16 | 155.34 | 155.37 | 155.41 | 155.52 | 155.52 | 155.53 | 155.56 | 155.57 | 155.57 | 155.86 | 155.86 | 155.87 | 156.13 | 156.21 | 156.51 | 156.51 | 156.90 | 157.03 |
| Delta | 0.00 | 0.02 | 0.07 | 0.12 | 0.31 | 0.34 | 0.38 | 0.49 | 0.49 | 0.50 | 0.53 | 0.53 | 0.54 | 0.83 | 0.83 | 0.83 | 1.10 | 1.18 | 1.47 | 1.48 | 1.87 | 1.99 |
| Weight | 0.03 | 0.03 | 0.03 | 0.03 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.02 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||||||||||||||||||||||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | Model 7 | Model 8 | Model 9 | |
|---|---|---|---|---|---|---|---|---|---|
| psi(Int) | 1.11 | 0.96 | 0.84 | 0.87 | 0.90 | 0.69 | 0.69 | 0.93 | 0.96 |
| (0.60) | (0.58) | (0.62) | (0.63) | (0.52) | (0.53) | (0.53) | (0.58) | (0.58) | |
| psi(Burn.intensity.basal) | 10.37 | 0.82 | 21.39 | 1.58 | |||||
| (8.19) | (0.83) | (18.95) | |||||||
| psi(Burn.intensity.Canopy) | -12.01 | 56.97 | -28.59 | 34.71 | 1.26 | ||||
| (9.25) | (208.37) | (24.56) | (92.31) | (1.94) | |||||
| p(Int) | 1.15 | 1.46* | 1.34* | 1.36* | 3.91 | 3.85 | 3.95 | 4.78 | 4.64 |
| (0.76) | (0.71) | (0.66) | (0.66) | (2.72) | (2.42) | (2.43) | (2.54) | (2.64) | |
| p(Meantemp) | -0.05 | -0.08 | -0.08 | -0.08 | -0.09 | -0.08 | |||
| (0.06) | (0.05) | (0.05) | (0.05) | (0.05) | (0.05) | ||||
| psi(Burn.intensity.soil) | 42.80 | 3.70 | 33.60 | ||||||
| (221.42) | (3.66) | (180.26) | |||||||
| p(Julian) | -0.02 | -0.02 | -0.02 | -0.02 | -0.02 | ||||
| (0.01) | (0.01) | (0.01) | (0.01) | (0.01) | |||||
| p(Meanhum) | 0.00 | ||||||||
| (0.01) | |||||||||
| Log Likelihood | -85.65 | -87.30 | -87.32 | -87.35 | -83.45 | -87.70 | -87.71 | -86.47 | -86.55 |
| AICc | 182.80 | 183.58 | 183.62 | 183.68 | 183.84 | 184.39 | 184.39 | 184.44 | 184.60 |
| Delta | 0.00 | 0.79 | 0.82 | 0.89 | 1.04 | 1.59 | 1.59 | 1.65 | 1.81 |
| Weight | 0.10 | 0.06 | 0.06 | 0.06 | 0.06 | 0.04 | 0.04 | 0.04 | 0.04 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | |||||||||
| Model 1 | Model 2 | |
|---|---|---|
| psi(Int) | 0.76 | 0.85 |
| (0.47) | (0.52) | |
| p(Int) | 5.86 | 5.03 |
| (3.00) | (3.04) | |
| p(Julian) | -0.03 | -0.03 |
| (0.02) | (0.02) | |
| p(Meanhum) | 0.01 | |
| (0.01) | ||
| Log Likelihood | -78.56 | -77.89 |
| AICc | 163.68 | 164.75 |
| Delta | 0.00 | 1.07 |
| Weight | 0.09 | 0.05 |
| Num. obs. | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | |
|---|---|---|---|---|---|---|
| psi(Int) | 8.31 | 7.98 | -0.45 | -0.08 | -0.36 | 5.89 |
| (89.85) | (65.83) | (2.28) | (1.17) | (2.35) | (21.66) | |
| p(Int) | 1.33 | -0.23 | 1.59 | 4.76 | 1.58 | 5.63 |
| (2.14) | (2.60) | (2.14) | (4.55) | (2.15) | (6.45) | |
| p(Meanhum) | -0.11 | -0.14 | -0.10 | -0.27 | -0.11 | -0.11 |
| (0.06) | (0.07) | (0.06) | (0.19) | (0.06) | (0.06) | |
| p(sdtemp) | 0.62 | |||||
| (0.50) | ||||||
| psi(Burn.intensity.basal) | 0.68 | |||||
| (0.90) | ||||||
| p(sdhum) | 0.35 | |||||
| (0.28) | ||||||
| psi(Burn.intensity.Canopy) | 0.87 | |||||
| (1.24) | ||||||
| p(Julian) | -0.02 | |||||
| (0.03) | ||||||
| Log Likelihood | -10.76 | -9.89 | -10.50 | -10.50 | -10.52 | -10.52 |
| AICc | 28.09 | 28.75 | 29.97 | 29.98 | 30.02 | 30.02 |
| Delta | 0.00 | 0.66 | 1.88 | 1.89 | 1.93 | 1.93 |
| Weight | 0.08 | 0.06 | 0.03 | 0.03 | 0.03 | 0.03 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||||||
| Model 1 | Model 2 | Model 3 | Model 4 | Model 5 | Model 6 | Model 7 | Model 8 | |
|---|---|---|---|---|---|---|---|---|
| psi(Int) | -1.83 | -1.25 | 3.25 | 7.75 | -1.55 | -0.94 | 2.31 | 3.93 |
| (1.69) | (1.25) | (3.34) | (10.64) | (1.50) | (1.57) | (3.90) | (4.07) | |
| psi(Burn.intensity.basal) | -137.65 | -90.37 | 73.50 | 216.34 | -139.14 | -82.45 | 108.19 | |
| (94.08) | (67.03) | (86.62) | (274.45) | (91.83) | (51.80) | (119.13) | ||
| psi(Burn.intensity.Canopy) | 137.00 | 90.09 | 138.17 | 83.54 | 95.68 | |||
| (92.80) | (65.42) | (90.24) | (51.62) | (113.69) | ||||
| p(Int) | -0.55 | 0.65 | -1.93*** | -0.80 | 0.64 | -1.88*** | -2.03*** | -2.50 |
| (1.75) | (1.53) | (0.45) | (1.15) | (1.62) | (0.46) | (0.44) | (1.40) | |
| p(Meantemp) | -0.23 | -0.17 | -0.21 | 0.06 | ||||
| (0.12) | (0.11) | (0.11) | (0.11) | |||||
| p(sdtemp) | 0.61 | |||||||
| (0.43) | ||||||||
| psi(Burn.intensity.soil) | -53.48 | -157.36 | -69.47 | -78.45 | ||||
| (63.22) | (199.55) | (83.89) | (86.50) | |||||
| p(Meanhum) | -0.01 | |||||||
| (0.02) | ||||||||
| p(sdhum) | 0.08 | |||||||
| (0.09) | ||||||||
| Log Likelihood | -15.21 | -16.90 | -18.27 | -17.03 | -15.99 | -18.63 | -18.67 | -17.53 |
| AICc | 44.58 | 45.31 | 45.45 | 45.56 | 46.13 | 46.17 | 46.24 | 46.56 |
| Delta | 0.00 | 0.73 | 0.87 | 0.98 | 1.55 | 1.59 | 1.66 | 1.98 |
| Weight | 0.09 | 0.06 | 0.06 | 0.05 | 0.04 | 0.04 | 0.04 | 0.03 |
| Num. obs. | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| p < 0.001, p < 0.01, p < 0.05 | ||||||||
## ~/Documents/new_bats/Rnew_bats/Bats_data_products/fire.asc has GDAL driver AAIGrid
## and has 250 rows and 434 columns
## ~/Documents/new_bats/Rnew_bats/Bats_data_products/not_fire.asc has GDAL driver AAIGrid
## and has 250 rows and 434 columns
| with.fire | without.fire | with.fire.sd | |
|---|---|---|---|
| Yuma.Myotis | 0.23 | 0.18 | 0.41 |
| Small.Footed.Myotis | 0.15 | 0.09 | 0.33 |
| Little.Brown.Bat | 0.13 | 0.21 | 0.26 |
| Western.Red.Bat | 0.12 | 0.53 | 0.31 |
| Long.eared.Bat | 0.22 | 0.65 | 0.39 |
| Pallid.Bat | 0.09 | 0.16 | 0.28 |
| Townsend.s.big.eared.Bat | 0.16 | 0.19 | 0.35 |
| Big.Brow.Bat | 0.21 | 0.35 | 0.39 |
| Silver.Haired.Bat | 0.18 | 0.49 | 0.36 |
| Brazilian.free.tailed.bat | 0.19 | 0.74 | 0.37 |
| western.mastiff.bat | 0.01 | 0.15 | 0.07 |
##
## One Sample t-test
##
## data: myyu
## t = 23.888, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.1818892
## 95 percent confidence interval:
## 0.2256631 0.2334887
## sample estimates:
## mean of x
## 0.2295759
##
## One Sample t-test
##
## data: myci
## t = 36.501, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.09200702
## 95 percent confidence interval:
## 0.1482991 0.1546876
## sample estimates:
## mean of x
## 0.1514933
##
## One Sample t-test
##
## data: mylu
## t = -57.863, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.2051132
## 95 percent confidence interval:
## 0.1304989 0.1353882
## sample estimates:
## mean of x
## 0.1329435
##
## One Sample t-test
##
## data: labl
## t = -264.62, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.5267009
## 95 percent confidence interval:
## 0.1164016 0.1224351
## sample estimates:
## mean of x
## 0.1194183
##
## One Sample t-test
##
## data: myev
## t = -229.73, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.6525984
## 95 percent confidence interval:
## 0.2128346 0.2202750
## sample estimates:
## mean of x
## 0.2165548
##
## One Sample t-test
##
## data: anpa
## t = -52.936, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.1632505
## 95 percent confidence interval:
## 0.08879395 0.09411075
## sample estimates:
## mean of x
## 0.09145235
##
## One Sample t-test
##
## data: coto
## t = -21.478, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.1948568
## 95 percent confidence interval:
## 0.1550577 0.1617141
## sample estimates:
## mean of x
## 0.1583859
##
## One Sample t-test
##
## data: epfu
## t = -75.202, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.3493452
## 95 percent confidence interval:
## 0.2018752 0.2093671
## sample estimates:
## mean of x
## 0.2056211
##
## One Sample t-test
##
## data: lano
## t = -173.63, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.4905386
## 95 percent confidence interval:
## 0.1775398 0.1845274
## sample estimates:
## mean of x
## 0.1810336
##
## One Sample t-test
##
## data: tabr
## t = -301.07, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.7393856
## 95 percent confidence interval:
## 0.1895036 0.1966171
## sample estimates:
## mean of x
## 0.1930604
##
## One Sample t-test
##
## data: eupe
## t = -413.34, df = 41875, p-value < 2.2e-16
## alternative hypothesis: true mean is not equal to 0.1519613
## 95 percent confidence interval:
## 0.008627646 0.009980571
## sample estimates:
## mean of x
## 0.009304108
library(vioplot)
## Loading required package: sm
## Package 'sm', version 2.2-5.4: type help(sm) for summary information
vioplot(myyu,myci, mylu, labl, myev, anpa, coto, epfu, lano, tabr, eupe, col="grey")
#valuetable <- getValues(AllLayers2)
#km1 <- kmeans(na.omit(valuetable), centers = 5, iter.max = 100, nstart = 10)
# create a blank raster with default values of 0
#rNA <- setValues(raster(AllLayers2), 0)
#for(i in 1:nlayers(AllLayers2)){
#rNA[is.na(AllLayers2[[i]])] <- 1
#}
# convert rNA to an integer vector
#rNA <- getValues(rNA)
# convert valuetable to a data.frame
#valuetable <- as.data.frame(valuetable)
# if rNA is a 0, assign the cluster value at that position
#valuetable$class[rNA==0] <- km1$cluster
# if rNA is a 1, assign an NA at that position
#valuetable$class[rNA==1] <- NA
# create a blank raster
#classes1 <- raster(AllLayers2)
# assign values from the 'class' column of valuetable
#classes1 <- setValues(classes1, valuetable$class)
#plot(classes1, legend=TRUE, colNA="black")
#More info on how to do this clasification in *https://geoscripting-wur.github.io/AdvancedRasterAnalysis/*
##
## Call:
## occu(formula = ~Meanhum + 1 ~ 1, data = SimOccuMyVo2)
##
## Occupancy:
## Estimate SE z P(>|z|)
## -0.33 0.731 -0.451 0.652
##
## Detection:
## Estimate SE z P(>|z|)
## (Intercept) -0.49767 1.1058 -0.450 0.653
## Meanhum -0.00923 0.0144 -0.641 0.522
##
## AIC: 85.66033
## 0.025 0.975
## p(Int) -2.66493507 1.66960349
## p(Meanhum) -0.03746837 0.01900186
## Profiling parameter 1 of 2 ... done.
## Profiling parameter 2 of 2 ... done.
## 0.025 0.975
## p(Int) -2.76769061 1.62413485
## p(Meanhum) -0.03895494 0.01837992
## Backtransformed linear combination(s) of Detection estimate(s)
##
## Estimate SE LinComb (Intercept) Meanhum
## 0.377 0.258 -0.502 1 0.5
##
## Transformation: logistic
## Length Class Mode
## 1 unmarkedBackTrans S4
## 0.05 0.95
## 0.09019551 0.786958
## Warning: Some observations have been discarded because corresponding
## covariates were missing.
## Warning: 8 sites have been discarded because of missing data.
## Warning: Some observations have been discarded because corresponding
## covariates were missing.
## Warning: 8 sites have been discarded because of missing data.
## [1] 10
##
## --------------------------------------------------------------------------
## Evaluation of design K = 4 S = 62 (TS = 248)
## --------------------------------------------------------------------------
## estimator performance (excl empty histories)
## psi: bias = +0.0310 var = +0.0148 MSE = +0.0158
## p: bias = -0.0087 var = +0.0102 MSE = +0.0103
## covar = -0.0075 critA = +0.0261 critD = +1.063e-04
## estimator performance (excl also histories leading to boundary estimates)
## psi: bias = +0.0219 var = +0.0079 MSE = +0.0084
## p: bias = -0.0055 var = +0.0095 MSE = +0.0095
## covar = -0.0051 critA = +0.0179 critD = +5.378e-05
## empty histories = 0.0%
## boundary estimates = 1.2%
## this took 1.576 seconds
## --------------------------------------------------------------------------
## $dist
## [,1] [,2] [,3] [,4] [,5]
## [1,] 11 18 0.0168 0.22452662 0.323314803
## [2,] 8 14 0.0155 0.15365326 0.367330550
## [3,] 6 7 0.0084 0.27681552 0.101952333
## [4,] 7 11 0.0235 0.14942362 0.296870257
## [5,] 6 13 0.0023 0.10264751 0.510559198
## [6,] 10 13 0.0121 0.30308670 0.172966754
## [7,] 13 20 0.0091 0.28505105 0.282991625
## [8,] 12 21 0.0097 0.23050924 0.367328088
## [9,] 6 11 0.0121 0.11138423 0.398138526
## [10,] 6 8 0.0142 0.17026449 0.189409474
## [11,] 17 26 0.0008 0.37564141 0.279046843
## [12,] 10 17 0.0185 0.19689844 0.348385401
## [13,] 11 19 0.0135 0.21341619 0.358850470
## [14,] 12 14 0.0019 0.55345996 0.102016630
## [15,] 8 10 0.0146 0.27368013 0.147316041
## [16,] 6 10 0.0150 0.12025816 0.335412443
## [17,] 7 9 0.0149 0.21899815 0.165661171
## [18,] 14 20 0.0050 0.34390571 0.234542948
## [19,] 11 17 0.0193 0.23980349 0.285788874
## [20,] 8 17 0.0035 0.13788019 0.497449536
## [21,] 13 16 0.0030 0.47047257 0.137147180
## [22,] 9 17 0.0111 0.16392186 0.418082881
## [23,] 12 16 0.0066 0.34057947 0.189403215
## [24,] 10 16 0.0218 0.20903637 0.308540833
## [25,] 15 20 0.0014 0.42589996 0.189426297
## [26,] 12 20 0.0115 0.24051946 0.335321201
## [27,] 4 6 0.0064 0.09085934 0.265426361
## [28,] 8 16 0.0067 0.14139524 0.456252348
## [29,] 16 24 0.0011 0.36327383 0.266357280
## [30,] 8 11 0.0217 0.21170026 0.209495206
## [31,] 9 14 0.0265 0.19452864 0.290223943
## [32,] 11 14 0.0078 0.35468138 0.159088653
## [33,] 12 13 0.0004 0.98663540 0.053153054
## [34,] 11 15 0.0126 0.29627610 0.204108286
## [35,] 9 16 0.0166 0.17077992 0.377800298
## [36,] 5 8 0.0114 0.10450134 0.308641291
## [37,] 6 9 0.0173 0.13620539 0.266398992
## [38,] 8 12 0.0271 0.18169610 0.266368078
## [39,] 8 13 0.0225 0.16437782 0.318776995
## [40,] 8 9 0.0049 0.46437881 0.078150915
## [41,] 14 24 0.0045 0.27340452 0.353871992
## [42,] 9 11 0.0099 0.33458757 0.132526116
## [43,] 13 19 0.0087 0.30712173 0.249400328
## [44,] 13 22 0.0081 0.25681489 0.345460150
## [45,] 4 5 0.0048 0.13692652 0.147275280
## [46,] 11 20 0.0082 0.20537835 0.392717080
## [47,] 11 16 0.0187 0.26187041 0.246388971
## [48,] 10 18 0.0124 0.18803852 0.386041268
## [49,] 9 18 0.0053 0.15890691 0.456408635
## [50,] 7 12 0.0178 0.13683439 0.353740348
## [51,] 12 18 0.0135 0.27242819 0.266428446
## [52,] 10 15 0.0230 0.22706994 0.266448455
## [53,] 6 12 0.0055 0.10603610 0.456406392
## [54,] 7 15 0.0025 0.12021707 0.503058970
## [55,] 5 9 0.0075 0.09400784 0.385927448
## [56,] 12 17 0.0123 0.29947289 0.228941738
## [57,] 11 21 0.0060 0.19922504 0.425185283
## [58,] 3 6 0.0019 0.05300865 0.456372065
## [59,] 5 7 0.0108 0.12757896 0.221400144
## [60,] 3 3 0.0014 1.00000000 0.012096774
## [61,] 15 26 0.0019 0.29027207 0.361130577
## [62,] 5 11 0.0017 0.08510886 0.521076113
## [63,] 8 15 0.0108 0.14645019 0.413210870
## [64,] 9 13 0.0233 0.21678472 0.241739377
## [65,] 15 19 0.0013 0.49172410 0.155752621
## [66,] 12 19 0.0147 0.25382091 0.301834517
## [67,] 9 12 0.0166 0.25534831 0.189575780
## [68,] 13 18 0.0064 0.33897054 0.214155235
## [69,] 13 21 0.0092 0.26899419 0.314836288
## [70,] 10 14 0.0189 0.25513422 0.221318134
## [71,] 2 3 0.0008 0.04543110 0.266305228
## [72,] 16 28 0.0016 0.30723307 0.367365899
## [73,] 6 14 0.0012 0.10052001 0.561841302
## [74,] 14 22 0.0065 0.29897145 0.296794096
## [75,] 12 24 0.0020 0.21207247 0.456215194
## [76,] 14 25 0.0033 0.26473166 0.380753640
## [77,] 9 21 0.0005 0.15075108 0.561834164
## [78,] 12 22 0.0057 0.22260995 0.398276496
## [79,] 9 15 0.0204 0.18039585 0.335309322
## [80,] 4 8 0.0029 0.07066999 0.456206586
## [81,] 10 22 0.0008 0.17026940 0.521078617
## [82,] 13 15 0.0008 0.63834760 0.094764169
## [83,] 16 23 0.0014 0.38869900 0.238559474
## [84,] 7 10 0.0201 0.17192187 0.234573551
## [85,] 15 23 0.0034 0.33035075 0.280689739
## [86,] 9 10 0.0027 0.57719723 0.069838114
## [87,] 14 23 0.0046 0.28455066 0.325980193
## [88,] 7 20 0.0001 0.11376986 0.709156825
## [89,] 10 19 0.0090 0.18155868 0.421805504
## [90,] 14 21 0.0054 0.31781398 0.266448182
## [91,] 7 13 0.0118 0.12887761 0.406942551
## [92,] 17 30 0.0004 0.32434912 0.372983805
## [93,] 12 15 0.0046 0.41056052 0.147318715
## [94,] 1 2 0.0001 0.01766963 0.456484625
## [95,] 4 7 0.0048 0.07679513 0.367546391
## [96,] 9 19 0.0015 0.15546292 0.492907481
## [97,] 7 14 0.0073 0.12373890 0.456442535
## [98,] 11 22 0.0037 0.19438207 0.456402507
## [99,] 7 8 0.0053 0.36480763 0.088405054
## [100,] 13 23 0.0058 0.24750599 0.374645093
## [101,] 14 18 0.0015 0.43805504 0.165678534
## [102,] 10 12 0.0061 0.40136511 0.120558967
## [103,] 14 26 0.0023 0.25773347 0.406709172
## [104,] 15 31 0.0001 0.26133983 0.478491389
## [105,] 17 28 0.0008 0.34471097 0.327540753
## [106,] 1 1 0.0003 1.00000000 0.004032258
## [107,] 15 25 0.0033 0.30061932 0.335393779
## [108,] 4 9 0.0013 0.06762562 0.536685098
## [109,] 17 23 0.0004 0.46573612 0.199047562
## [110,] 15 29 0.0005 0.26957139 0.433753934
## [111,] 5 10 0.0038 0.08831069 0.456447998
## [112,] 5 6 0.0077 0.20071554 0.120553856
## [113,] 7 7 0.0014 1.00000000 0.028225806
## [114,] 14 27 0.0009 0.25207379 0.431944947
## [115,] 10 21 0.0018 0.17310340 0.489276407
## [116,] 3 4 0.0023 0.08510598 0.189577897
## [117,] 12 23 0.0035 0.21682526 0.427818610
## [118,] 16 22 0.0006 0.42343824 0.209490427
## [119,] 12 25 0.0011 0.20836008 0.483787000
## [120,] 10 11 0.0015 0.70222585 0.063152222
## [121,] 4 4 0.0021 1.00000000 0.016129032
## [122,] 5 5 0.0010 1.00000000 0.020161290
## [123,] 17 32 0.0003 0.31036046 0.415846762
## [124,] 11 23 0.0010 0.19066611 0.486364326
## [125,] 3 8 0.0001 0.04908538 0.657482614
## [126,] 17 25 0.0005 0.39769141 0.253463818
## [127,] 13 17 0.0050 0.38773399 0.176805026
## [128,] 9 20 0.0007 0.15273405 0.527874018
## [129,] 15 22 0.0024 0.35245447 0.251647544
## [130,] 14 19 0.0020 0.38100614 0.201063680
## [131,] 11 13 0.0031 0.47452880 0.110444190
## [132,] 16 26 0.0018 0.32895778 0.318760744
## [133,] 9 9 0.0007 1.00000000 0.036290323
## [134,] 3 7 0.0003 0.05022243 0.561871836
## [135,] 20 27 0.0001 0.55128443 0.197552156
## [136,] 9 22 0.0001 0.14919924 0.594665446
## [137,] 11 24 0.0006 0.18768233 0.515639224
## [138,] 18 26 0.0003 0.43362577 0.241766152
## [139,] 7 16 0.0010 0.11781088 0.547414033
## [140,] 6 6 0.0025 1.00000000 0.024193548
## [141,] 8 8 0.0013 1.00000000 0.032258065
## [142,] 13 26 0.0013 0.22972210 0.456282298
## [143,] 15 30 0.0002 0.26505976 0.456327633
## [144,] 14 17 0.0007 0.53430306 0.128314902
## [145,] 17 27 0.0007 0.35850592 0.303687196
## [146,] 3 5 0.0013 0.06009815 0.335485005
## [147,] 15 27 0.0014 0.28190827 0.385994888
## [148,] 18 31 0.0002 0.35021189 0.356964028
## [149,] 18 29 0.0002 0.37351872 0.313087763
## [150,] 16 25 0.0017 0.34398051 0.293100674
## [151,] 15 24 0.0036 0.31364246 0.308595169
## [152,] 8 18 0.0015 0.13525324 0.536629536
## [153,] 17 22 0.0001 0.52194785 0.169925150
## [154,] 13 24 0.0027 0.24011773 0.402867929
## [155,] 10 20 0.0027 0.17669887 0.456345078
## [156,] 2 4 0.0003 0.03532969 0.456304219
## [157,] 16 27 0.0008 0.31695999 0.343549308
## [158,] 15 21 0.0012 0.38251227 0.221464229
## [159,] 16 29 0.0004 0.29937922 0.390620531
## [160,] 17 24 0.0007 0.42679485 0.226759005
## [161,] 14 28 0.0002 0.24741321 0.456344268
## [162,] 13 25 0.0023 0.23437686 0.430032162
## [163,] 15 18 0.0004 0.60239444 0.120436207
## [164,] 15 28 0.0006 0.27519917 0.410221285
## [165,] 11 12 0.0006 0.83791891 0.057766725
## [166,] 2 2 0.0005 1.00000000 0.008064516
## [167,] 5 13 0.0003 0.08203616 0.639079748
## [168,] 19 27 0.0001 0.47125747 0.231027034
## [169,] 16 21 0.0008 0.47282254 0.179137799
## [170,] 6 15 0.0004 0.09910886 0.610637684
## [171,] 11 11 0.0001 1.00000000 0.044354839
## [172,] 2 5 0.0002 0.03302487 0.610640212
## [173,] 12 29 0.0001 0.19936828 0.586559302
## [174,] 18 30 0.0006 0.36066953 0.335321309
## [175,] 13 27 0.0008 0.22595535 0.481784272
## [176,] 18 32 0.0001 0.34153464 0.377888136
## [177,] 13 14 0.0001 1.00000000 0.056451613
## [178,] 14 16 0.0002 0.72911735 0.088482574
## [179,] 7 18 0.0001 0.11503092 0.630893043
## [180,] 12 26 0.0003 0.20530490 0.510634339
## [181,] 17 29 0.0006 0.33344069 0.350621546
## [182,] 19 29 0.0001 0.42108035 0.277624695
## [183,] 12 27 0.0002 0.20289593 0.536513755
## [184,] 14 30 0.0002 0.24041393 0.503033379
## [185,] 16 30 0.0004 0.29282014 0.413250184
## [186,] 20 30 0.0001 0.45408532 0.266343361
## [187,] 17 21 0.0001 0.60716623 0.139481193
## [188,] 10 10 0.0002 1.00000000 0.040322581
## [189,] 6 16 0.0001 0.09813565 0.657551647
## [190,] 13 29 0.0002 0.22037087 0.530679657
## [191,] 5 12 0.0005 0.08319267 0.581545155
## [192,] 19 26 0.0001 0.50790235 0.206386984
## [193,] 10 23 0.0001 0.16804222 0.551718006
## [194,] 18 27 0.0002 0.40884225 0.266290847
## [195,] 17 33 0.0001 0.30498893 0.436340091
## [196,] 18 25 0.0001 0.46639830 0.216092923
## [197,] 7 17 0.0002 0.11619216 0.589896956
## [198,] 18 28 0.0004 0.38924195 0.290088806
## [199,] 20 34 0.0001 0.39354030 0.348300999
## [200,] 19 32 0.0001 0.37695198 0.342237148
## [201,] 4 11 0.0001 0.06521474 0.680227931
## [202,] 4 10 0.0004 0.06599048 0.610636843
## [203,] 8 19 0.0001 0.13337892 0.574232112
## [204,] 20 35 0.0001 0.38406358 0.367423642
## [205,] 22 32 0.0001 0.52398447 0.246280754
## [206,] 16 34 0.0001 0.27571023 0.497370046
## [207,] 18 36 0.0001 0.31821365 0.456263283
## [208,] 18 23 0.0001 0.57428424 0.161588753
## [209,] 21 34 0.0001 0.43342594 0.316365050
## [210,] 19 34 0.0001 0.35862602 0.382214904
## [211,] 10 24 0.0001 0.16640443 0.581549235
## [212,] 16 20 0.0001 0.54739008 0.147306892
## [213,] 19 31 0.0001 0.38898337 0.321408794
## [214,] 13 28 0.0001 0.22286876 0.506517918
## [215,] 19 30 0.0001 0.40340729 0.299863283
## [216,] 14 29 0.0001 0.24359192 0.480005463
## [217,] 12 12 0.0001 1.00000000 0.048387097
## [218,] 18 24 0.0001 0.51085751 0.189456468
##
## $biaspsi
## [1] 0.03102304
##
## $varpsi
## [1] 0.014844
##
## $MSEpsi
## [1] 0.01580643
##
## $biasp
## [1] -0.008716855
##
## $varp
## [1] 0.01019557
##
## $MSEp
## [1] 0.01027155
##
## $covar
## [1] -0.007486403
##
## $critA
## [1] 0.02607799
##
## $critD
## [1] 0.0001063104
##
## $biaspsi_B
## [1] 0.02191949
##
## $varpsi_B
## [1] 0.007936446
##
## $MSEpsi_B
## [1] 0.00841691
##
## $biasp_B
## [1] -0.005538184
##
## $varp_B
## [1] 0.009451622
##
## $MSEp_B
## [1] 0.009482293
##
## $covar_B
## [1] -0.005101769
##
## $critA_B
## [1] 0.0178992
##
## $critD_B
## [1] 5.378356e-05
##
## $pempty
## [1] 0
##
## $pbound
## [1] 1.17